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Abstract 


The first search for supersymmetry in the vector-boson fusion topology is presented. 
The search targets final states with at least two leptons, large missing transverse 
momentum, and two jets with a large separation in rapidity The data sample 
corresponds to an integrated luminosity of 19.7 fb^^ of proton-proton collisions at 
\/s = 8 TeV collected with the CMS detector at the CERN LHC. The observed dijet 
invariant mass spectrum is found to be consistent with the expected standard model 
prediction. Upper limits are set on the cross sections for chargino and neutralino pro¬ 
duction with two associated jets, assuming the supersymmetric partner of the t lep¬ 
ton to be the lightest slepton and the lightest slepton to be lighter than the charginos. 
For a so-called compressed-mass-spectrum scenario in which the mass difference be¬ 
tween the lightest supersymmetric particle ^ and the next lightest, mass-degenerate, 
gaugino particles ^ and Xi" is 50 GeV, a mass lower limit of 170 GeV is set for these 
latter two particles. 
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1 Introduction 

With the successful operation of the CERN LHC, numerous results placing constraints on ex¬ 
tensions to the standard model (SM) have been presented by the ATLAS and CMS experi¬ 
ments. In particular, in models of supersymmetry (SUSY) lHHZl, limits in excess of 1 TeV have 
been placed on the masses of the strongly produced gluinos and first- and second-generation 
squarks UStiTSt. In contrast, mass limits on the weakly produced charginos (xf) and neutrali- 
nos i^), with much smaller production cross sections, are much less severe. The limits for 
charginos and neutralinos are especially weak in so-called compressed-mass-spectrum scenar¬ 
ios, in which the mass of the lightest supersymmetric particle (LSP) is only slightly less than 
the masses of other SUSY states. The chargmo-neutralmo sector plays a crucial role in the con¬ 
nection between dark matter and SUSY: in SUSY models with R-parity |T^ conservation, the 
lightest neutralino ^ often takes the role of the LSP and is a dark matter candidate. 

Previous LHC searches IITtI [TSl for electroweak chargino and neutralino production have fo¬ 
cused on final states with one or more leptons {£) and missing transverse momentum 
e.g., pair production followed by x^ —> £v^ and ^ > ££^, where x^ (^) is the light¬ 

est (next-to-lightest) chargino (neutralino), and where the LSP ^ is presumed to escape with¬ 
out detection leading to significant However, these searches exhibit limited sensitivity 

in cases where the x^ and ^ are nearly mass degenerate with the The mass difference 
Am = m~± ~ is a crucial parameter dictating the sensitivity of the analysis. While the ex¬ 
clusion limits in Refs. lIZllISl can be as large as m~± < 720 GeV for a massless they weaken 

to only « 100 GeV for Am < 50 GeV. The current searches also exhibit limited sensitivity to 
models with SUSY particles that decay predominantly to t leptons, even for LSP masses near 
zero, due to the larger backgrounds associated with r-lepton reconstruction compared to elec¬ 
trons or muons. 

Electroweak SUSY particles can be pair produced in association with two jets in pure elec¬ 
troweak processes in the vector-boson fusion (VBL) topology 1191, which is characterized by 
the presence of two forward jets (i.e. jets near the beam axis), in opposite hemispheres, lead¬ 
ing to a large dijet invariant mass {mjj). Ligure shows the Leynman diagrams for two of the 
possible VBL production processes: chargino-neutralino and chargino-chargino production. 



Ligure 1: Diagrams of (left) chargino-neutralino and (right) chargino-chargino pair production 
through vector-boson fusion followed by their decays to leptons and the LSP 


A search in the VBL topology offers a new and complementary means to directly probe the 
electroweak sector of SUSY, especially in compressed-mass-spectrum scenarios Il20l . It targets 
unexplored regions of SUSY parameter space, where other searches have limited sensitivity. 
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2 CMS detector 


It differs fundamentally from the conventional direct electroweak SUSY searches mentioned 
above in that it uses the presence of jets with large transverse momenta (px) to suppress SM 
background. In this regard, it resembles searches for strongly produced SUSY particles. How¬ 
ever, unlike these latter studies, which present searches for the indirect production of charginos 
and neutralinos through squark or gluino decay chains lUOl - tT^ . the VBF search does not require 
the production of squarks or gluinos, whose masses might be too large to allow production at 
the LHC. 

In this paper, we present a search for the electroweak production of SUSY particles in the VBF 
topology. The data, corresponding to an integrated luminosity of 19.7fb^^ of proton-proton 
collisions at a centre-of-mass energy of = 8 TeV, were collected with the CMS detector in 
2012. Besides the two oppositely directed forward jets (/) that define the VBF configuration, the 
search requires the presence of at least two leptons (e, p, or r) and large The events are 

classified into one of eight final states depending on the dilepton content and charges ep;;, pp/// 
pTh/;, and TtiThjj, where Th denotes a hadronically decaying t lepton and where we differen¬ 
tiate between final states with same-sign (SS) and opposite-sign (OS) dilepton pairs. The dijet 
invariant mass distribution nijj is used to search for the SUSY signal. Stringent requirements are 
placed on p^'®® and on the kinematic properties of the VBF dijet system to suppress SM back¬ 
ground. In particular, the R-parity conserving SUSY models we examine predict much higher 
average dijet px than is typical for SM processes such as VBF Higgs boson production EH . 
allowing us to suppress the background by a factor of 10^-10^, depending on the background 
process. 

The background is evaluated using data wherever possible. The general strategy is to define 
control regions, each dominated by a different background process, through modification of 
the nominal selection requirements. These control regions are used to measure the nijj shapes 
and probabilities for background events to satisfy the VBF selection requirements. If the back¬ 
ground contribution from a particular process is expected to be small or if the above approach 
is not feasible, the nijj shapes are taken from simulation. In these cases, scale factors, defined as 
the ratio of efficiencies measured in data and simulation, are used to normalize the predicted 
rates to the data. 

The paper is organized as follows. The CMS detector is described in Section!^ The reconstruc¬ 
tion of electrons, muons, Th leptons, jets, and pi^^®® is presented in Section^ The dominant 
backgrounds and their simulated samples are discussed in Section followed by the descrip¬ 
tion of the event selection in Section]^ and the background estimation in Section]^ Systematic 
uncertainties are summarized in Section]^ and the results are presented in Section]^ Section 9 
contains a summary. 


2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diame¬ 
ter, providing a magnetic field of 3.8 T. Located within the solenoid volume are a silicon pixel 
and strip tracker, a lead tungstate electromagnetic calorimeter (ECAL), and a brass and scin¬ 
tillator hadron calorimeter. Muons are measured in gas-ionisation detectors embedded in the 
steel flux-return yoke outside the solenoid. Extensive forward calorimetry complements the 
coverage provided by the barrel and endcap detectors. Forward hadron calorimeters on each 
side of the CMS interaction point cover the very forward angles of CMS, in the pseudorapidity 
range 3.0 < IpI < 5.0. 

The irmer tracker measures charged tracks with |p| < 2.5 and provides an impact parame- 
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ter resolution of ftslSf/m and a transverse momentum resolution of about 1.5% for 100 GeV 
charged particles. Events are selected with a first-level trigger made of a system of fast elec¬ 
tronics, and a high-level trigger that consists of a farm of commercial CPUs running a version 
of the offline reconstruction optimized for fast processing. A detailed description of the CMS 
detector, along with a definition of the coordinate system and relevant kinematic variables, can 
be found in Ref. I S21 . 


3 Object reconstruction and identification 

The missing transverse momentum vector is defined as the projection on the plane per¬ 
pendicular to the beam axis of the negative vector sum of the momenta of all reconstructed 
particles in an event. Its magnitude is referred to as p^^®®. The jets and p!p^®® are reconstructed 
with the particle-flow algorithm Il23ll^ . The anti-kj clustering algorithm f25l with a distance 
parameter of 0.5 is used for jet clustering. Jets are required to satisfy identification criteria 
designed to reject particles from multiple proton-proton interactions (pileup) and anomalous 
behavior from the calorimeters. For jets with pj > 30GeV and \rj\ < 2.5 (2.5 < \n\ < 5.0), 
the reconstruction-plus-identification efficiency is «99% (95%), while 90-95% (60%) of pileup 
jets are rejected 11261 . The jet energy scale and resolution are calibrated through correction fac¬ 
tors that depend on the px and p of the jet [27 1. Jets originating from the hadronisation of 
bottom quarks (b quark jets) are identified using the loose working point of the combined sec¬ 
ondary vertex (CSV) algorithm (281, which exploits observables related to the long lifetime of 
b hadrons. For jets with px > 20 GeV and |p| < 2.4, the probability of correctly identifying a b 
quark jet is ~85%, while the probability of misidentifying a jet originating from a light quark 
or gluon as a b quark jet is ?al0% [291 . 

Muons are reconstructed (301 using the inner silicon tracker and muon detectors. Quality re¬ 
quirements based on the minimum number of hits in the silicon tracker, pixel detector, and 
muon detectors are applied to suppress backgrounds from decays-in-flight and hadron shower 
remnants that reach the muon system. Electrons are reconstructed (311 by combining tracks 
produced by the Gaussian-sum filter algorithm with ECAF clusters. Requirements on the 
track quality, the shape of the energy deposits in the ECAL, and the compatibility of the mea¬ 
surements from the tracker and the ECAL are imposed to distinguish prompt electrons from 
charged pions and from electrons produced by photon conversions. The electron and muon 
reconstruction efficiencies are >99% for px > 10 GeV. 

The electron and muon candidates are required to satisfy isolation criteria in order to reject 
non-prompt leptons from the hadronisation process. Isolation is defined as the scalar sum of 
the Px values of reconstructed charged and neutral particles within a cone of radius AR = 
\/ {Apf- -h (A^)^ = 0.4 around the lepton-candidate track, divided by the px of the lepton can¬ 
didate. A correction is applied to the isolation variable to account for the effects of additional 
interactions. For charged particles, only tracks associated with the primary vertex are included 
in the isolation sums. The primary vertex is the reconstructed vertex with the largest sum of 
charged-track pj values associated to it. For neutral particles, a correction is applied by sub¬ 
tracting the energy deposited in the isolation cone by charged particles not associated with the 
primary vertex, multiplied by a factor of 0.5. This factor corresponds approximately to the ratio 
of neutral to charged hadron production in the hadronisation process of pileup interactions. In 
both cases, the contribution from the electron or muon candidate is removed from the sum and 
the value of the isolation variable is required to be 0.2 or less. 

The muon identification-plus-isolation efficiency is 96% for muons with Px > 15 GeV and |p| < 
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4 Signal and background samples 


2.1. The rate at which pions undergoing tt^ —> /r decay are misidentified as muons is 10^^ for 
pions with px > lOGeV and |p| < 2.1. The electron identification-plus-isolation efficiency 
is 85% (80%) for electrons with pj > 30GeV in the barrel (endcap) region lH]. The ; —> e 
misidentification rate is 5 x 10^^ for jets with px > lOGeV and |p| < 2.1. 

Hadronic decays of t leptons are reconstructed and identified using the hadrons-plus-strips 
algorithm 11321 , which is designed to optimize the performance of the Th reconstruction by con¬ 
sidering specific Th decay modes. To suppress backgrounds in which light-quark or gluon jets 
can mimic decays, a candidate is required to be spatially isolated from other energy de¬ 
posits in the event. The isolation variable is calculated using a multivariate boosted-decision- 
tree technique by forming rings of radius AR around the direction of the candidate, using 
the energy deposits of particles not considered in the reconstruction of the Th decay mode. Ad¬ 
ditionally, Th candidates are required to be distinguishable from electrons and muons in the 
event by using dedicated criteria based on the consistency among the measurements in the 
tracker, calorimeters, and muon detectors. The identification and isolation efficiency is 55-65% 
for a Th lepton with px > 20GeV and |p| < 2.1, depending on the px and rj values of the Th 
candidate. The rate at which jets are misidentified as a Th lepton is 1-5%, depending on the px 
and p values of the Th candidate. 

The event selection criteria used in each search channel are summarized in Section (see Ta¬ 
ble [^. 


4 Signal and background samples 

The composition of SM background events depends on the final state and, in particular, the 
number of Th candidates. The most important sources of background arise from the production 
of W or Z bosons in association with jets (W/ Z+jets), and from tt, diboson (VV: WW, WZ, ZZ), 
and Quantum GhromoDynamics (QGD) multijet events. The W+jets events are characterized 
by an isolated lepton from the decay of the W boson and uncorrelated jets misidentified as 
an e, p, or Th. Background from W+jets events is especially pertinent for final states with one 
Th candidate. Background from tf events usually contains one or two tagged b quark jets, in 
addition to a genuine isolated e, p, or Th. 

Background from diboson events contains genuine, isolated leptons when the bosons decay 
leptonically, and jets that are misidentified as a Th lepton when they decay hadronically. The 
QGD background is characterized by jets that are misidentified as an e, p, or Th lepton. The 
QGD multijet process is an appreciable background only for the ThTh final states. 

There are negligible contributions from background processes such as single-top and VBF pro¬ 
duction of a Higgs or Z boson. These background yields are taken from simulation. 

Simulated samples of signal and background events are generated using Monte Garlo (MG) 
event generators. The signal event samples are generated with the MadGraph program 
(v5.1.5) |33l, considering pair production of xf and ^ gauginos (XiXh and ^^) 

with two associated partons. The signal events are generated requiring a pseudorapidity gap 
|Ap| > 4.2 between the two partons, with px > 30 GeV for each parton. Background event sam¬ 
ples with a Higgs boson produced through VBF processes, and single top are generated with 
the POWHEG program (vl.0rl380) l34l . The MadGraph generator (v5.1.3) is used to describe 
Z+jets, W+jets, tf, diboson, and VBF Z boson production. 

The MG background and signal yields are normalized to the integrated luminosity of the data. 
The tf background is normalized to the next-to-next-to-leading-logarithm level using the cal- 
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culations of Refs. Il35ll^ . The Z+jets and W+jets processes are normalized to next-to-next-to- 
leading-order using the results from the FEWZ v2.1 Il37ll generator. The diboson background 
processes are normalized to next-to-leading-order using the MCFM v5.8 l38l generator, while 
the VBF Z boson events are normalized to next-to-leadmg order using the VBFNLO (v2.6) II39P 
] program. The single-top and VBF Higgs boson background yields are taken from the powheg 
program, where the next-to-leading order effects are incorporated. Signal cross sections are 
calculated at leading order using the MadGraph generator. All MC samples incorporate the 
CTEQ6L1 I40l or CTEQ6M f41] parton distribution functions (PDF). The corresponding evalu¬ 
ation of uncertainties in the signal cross sections is discussed in Section The range of signal 
cross sections is 50-1 fb for ^ masses of 100-400 GeV. The POWHEG and MadGraph 

generators are interfaced with the pythia (v6.4.22) 1421 program, which is used to describe the 
parton shower and hadronisation processes. The decays of t leptons are described using the 
TAUOLA (27.1215) Il43l program. 

The background samples are processed with a detailed simulation of the GMS apparatus using 
the Geant4 package [Ml/ while the response for signal samples is modeled with the GMS fast 
simulation program Il45l . For the signal acceptance and myy shapes based on the fast simulation, 
the differences with respect to the GEANT4-based results are found to be small (<5%). Gorrec- 
tions are applied to account for the differences. For all MG samples, multiple proton-proton in¬ 
teractions are superimposed on the primary collision process, and events are reweighted such 
that the distribution of reconstructed collision vertices matches that in data. The distribution 
of the number of pileup interactions per event has a mean of 21 and a root-mean-square of 5.5. 


5 Event selection 

A single-muon trigger Il30l with a pj threshold of 24 GeV is used for the ppjj, epjj, and }iT\^jj 
final states. The ThT^jj channels use a double-Th trigger 14^ that requires pj > 35 GeV for each 
Th. A requirement on pseudorapidity (\ri\ < 2.1) is applied to select high quality and well- 
isolated leptons (e, p, Th) within the tracker acceptance. The pj thresholds defining the search 
regions are chosen to achieve a trigger efficiency greater than 90%. For final states with at least 
one muon ippjj, epjj, pThjj), events are selected by requiring a muon with pj > 30 GeV. For 
the ThTh channels, both Th candidates are required to satisfy pj > 45 GeV. 

The following requirements are referred to as the "central selection", and are applied to all 
final states. Pairs of leptons are required to be separated by AR > 0.3 and to originate from the 
primary vertex. All charmels require exactly two leptons satisfying selection criteria. Events 
with an e or ^ are required to have > 75 GeV, while a requirement > 30 GeV is used 
for the ThTh/y final state to compensate for the loss in acceptance due to the higher pj threshold 
of Th leptons while maintaining similar background rejection. Background from tf events is 
reduced by removing events in which any jet has pj > 20 GeV, is separated from the leptons by 
AR > 0.3, and is identified as b-quark jet using the loose working point of the GSV algorithm. 

The "VBF selection" refers to the requirement of at least two jets in opposite hemispheres 
ipipi < 0) with large separation (|Ap| > 4.2). Events are selected with at least two jets with 
Pj > 50 GeV and pseudorapidity \p\ < 5.0. The search region has a lower background 

rate with respect to other final states, which makes it possible to relax the jet pj requirement 
to 30 GeV to increase the signal acceptance. The event selection criteria with pj > 30 GeV are 
referred to as "Loose". The event selection criteria with pj > 50 GeV are referred to as "Tight". 
Selected events are required to have a dijet candidate with nijj > 250 GeV. 

The signal region (SR) is defined as the events that satisfy the central and VBF selection criteria. 
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6 Background estimation 


A summary of the event selection criteria used in each channel is presented in Table 

Table 1: Summary of the event selection criteria for the different final states. The selections 
for the and e}ijj channels are presented in one column i£e/fiHjj) they are similar. The 
symbol ^e,f/,Th means that the lepton could be an electron, a muon, or a lepton. 


Selection 



ThTh/; 

pj{}i)[GeV] 

>30 

>30 

— 

pT(4/,)[GeV] 

>15(e), >10(f/) 

— 

— 

PT(Th)[GeV] 

— 

>20 

>45 

^ (4+'Th ) 

<2.1 

<2.1 

<2.1 

jets 

0 

0 

0 

p^i®®[GeV] 

>75 

>75 

>30 

Ft (jets) 

>30/50 

>50 

>30 

If (jets)! 

<5 

<5 

<5 

|A? 7 (jets)| 

>4.2 

>4.2 

>4.2 

FiF2 

<0 

<0 

<0 


6 Background estimation 

The general methodology used to evaluate the background is the same for all final states. We 
isolate various control regions (CR) to measure the VBF efficiencies (the probability for a given 
background component to satisfy the VBF selection criteria) and nijj shapes from data, validate 
the modeling of the central selection criteria, and determine a correction factor to account for 
the selection efficiency by assessing the level of agreement between data and simulation. For 
each final state, the same trigger is used for the CRs as for the corresponding SR. The VBF 
efficiency, measured in a CR satisfying only the central selection, is defined as the fraction of 
events in the CR additionally passing the VBF event selection criteria. 

The tt, W+jets, and VV backgrounds are evaluated using the following equation: 

= N^^(central) SF^|\central) (1) 

where N^q‘^ is the predicted background yield in the signal region, (central) is the rate 
predicted by the "BG" simulation (with BG = tt, W+jets, or VV) for the central selection, 
(central) is the data-to-simulation correction factor for the central region, given by the 
ratio of data to the "BG" simulation in control region GRl, and is the VBF efficiency, de¬ 
termined as a function of rrijj in data control sample GR2 or, in the case of VV events, from 
simulation. 

The event selection criteria used to define the GR must not bias the rrijj distribution. This is 
verified, in simulation and data, by comparing the rrijj distributions with and without the se¬ 
lection criteria used to define the GR. The background estimation technique used to measure 
the VBF efficiency and rrijj shape from data is performed with simulated events to test the clo¬ 
sure, where closure refers to the ability of the method to predict the correct background yields 
when using simulation in place of data. The closure tests demonstrate that the background 
determination techniques, described in detail below, reproduce the expected background dis¬ 
tributions in both rate and shape to within the statistical uncertainties. The difference between 
the nominal MG background yields and the yields predicted from the closure test are added 
in quadrature with the statistical uncertainties of the prediction to define a systematic uncer¬ 
tainty. Simulated samples are further used to verify that the composition of objects erroneously 
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identified as leptons, and their kinematic properties, are similar between the CRs and SR, and 
thus that the scale factors (central) are unbiased. A variety of generators (MadGraph, 
PYTHIA, and powheg) are used for this purpose to establish the robustness of this expectation. 

The production of tf events is an important source of background for the mijj, e}ijj, and }iT\^jj 
final states. Control regions enriched with tf events are obtained by requiring the presence of at 
least one reconstructed b-tagged jet with pj > 20 GeV. The tf purity of the resulting data CRls 
depends on the final state, ranging from 76 to 99%. The contributions of backgrounds other 
than tf events are subtracted from the data CRls using simulation. The scale factors are 

then determined. The uncertainties related to the subtraction procedure are propagated to the 
scale factors. For the OS channels, the scale factors are consistent with unity to within 3%. The 
tf events with OS lepton pairs arise from genuine isolated leptons produced through leptonic 
W boson decay, and are well modeled by the simulation. On the other hand, tf events with 
SS lepton pairs mostly contain a lepton candidate that is a misidentified hadron or jet, which 
is more difficult to accurately simulate. The scale factors for SS events range from 1.2 to 1.5, 
with statistical uncertainties up to 25%. Since the fraction of lepton candidates that are in fact a 
misidentified hadron or jet varies according to the final state, the scale factors are determined 
independently for each channel. In contrast, the VBF efficiency for a given combination of lep¬ 
ton flavors does not depend on the charge state, and thus each pair of final states with the same 
flavor combination shares the same VBF efficiency value. The VBF efficiency is measured in 
data GR2 control samples obtained by additionally removing the charge requirement on the 
lepton pair and relaxing or inverting the lepton isolation requirement (isolation sidebands) in 
order to enhance the purity and statistical precision of the GR2s. Figure shows the "Tight" 
and "Loose" VBF efficiencies measured from data, as a function of nijj, for events in the tf 
GR2s of the ppjj final states. The VBF efficiencies range from 0.02 to 0.003, with relative 
uncertainties below 11% for nijj > 250 GeV. We verify that the b tagging, charge, and isola¬ 
tion requirements used to obtain the GRl and GR2 samples do not bias the nijj shape or the 
kinematic distributions of the leptons. 

The production of W+jets events presents an important source of background only for the pThjj 
search channels. Samples enriched in W+jets events, with about 70% purity according to sim¬ 
ulation, are obtained by requiring the transverse mass ntj = \/2p!JT®®pj^[l — cos(A^^j ^miss )] be¬ 
tween and the muon transverse momentum pj to satisfy 40 < nij <110 GeV. The correc¬ 
tion factor SF^^-gjg is determined to be 0.90 ± 0.11, where the uncertainty is a combination of 
the statistical uncertainty from data, the statistical uncertainty from simulation, and the system¬ 
atic uncertainty associated with the subtraction of the non-W+jets backgrounds from the data 
control sample. The lepton and Th isolation sidebands are used to obtain W+jets-enriched GR2 
samples, with negligible expected signal contributions, not only to measure the VBF efficien¬ 
cies and nijj shape from data, but also to provide further validation of the correction 

factor. To validate the correction factor, the W+jets rate in the isolation sideband is scaled 
by 0.90 ± 0.11, and agreement between the data and the corrected W+jets prediction from sim¬ 
ulation is observed. The VBF efficiency determined from the GR2 control sample is measured 
to be ?al% for nijj > 250 GeV. Agreement between the VBF efficiencies of Z+jets and W+jets 
processes is observed in the channel. 

The background from VV events is significant for final states containing only electrons and 
muons, comprising up to 10% of the total SM background in the OS channels, and up to 40% 
in the SS charmels. The diboson background is suppressed in the Tk final states because of 
the lower average pj of the visible r-lepton decay products. Diboson events have genuine 
isolated leptons and p!p^®® from neutrinos and can satisfy the VBF selection when produced in 
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6 Background estimation 



Figure 2: The VBF efficiency (see text) as a function of jet pair mass nijj measured for the tf 
and Z+jets control regions of the mijj final state, for the "Loose" {pj > 30GeV) and "Tight" 
(Pt > SOGeV) event selections. 

association with jets arising from initial-state radiation or from a SM VBF process. We select 
diboson-enriched GRl samples (97% purity) by requiring at least three leptons and inverting 
the requirement (p“'®® < 75GeV). The level of agreement between data and simulation 
for the event rates, VBF efficiencies, and nijj shapes are found to be the same for all types of W 
processes in the GRl samples. The data-to-simulation correction factor is = 1.12 ± 0.06. 

The nijj distributions, following the VBF selections, are consistent between data and simulation. 
Therefore, the VBF efficiency is taken directly from simulation. 

The Z+jets background is important for all OS final states. This background is evaluated using 
the following relation: 

= ^&-fets(central) (central) evBk^jj)^ (2) 

where and have the same meaning as the analogous quantities in Eq. (jl 

(with BG = Z+jets), and the term is described below. Gontrol samples (GRl) dominatec 

by Z —> ^£+jets events with £ = e, p (> 98% purity) are selected by requiring p™‘®® < 75GeV 
and an OS lepton pair mass m^i consistent with the Z boson (60 < < 120 GeV). The rates 

and kinematic distributions of leptons in these control samples are consistent between the data 
and simulation. These control samples are used to determine both the SF^F^^^ correction fac¬ 
tors and the Cyf p terms, in the same manner as described above for the analogous quantities. 
The correction factors are unity to within 1%. Figure 2 shows the "Tight" and "Loose" VBF 
efficiencies measured from data, as a function of nijj, for events in the ppjj channel. The VBF 
efficiencies range from 10^^ to 10^^. The measured Cyfp terms agree with the results from 
simulation within 23%, which is taken as a systematic uncertainty both for the background 
prediction and for the VBF efficiency in simulated signal events. Additional orthogonal Z+jets 
control samples (GR3) are selected with similar selection criteria as used for the signal, main¬ 
taining the p^^®® requirement (>75GeV) and inverting the VBF selection (i.e. at least one of the 
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VBF selection requirements is not satisfied: >2 jets, jet pj, \Ar]\, or riit] 2 )- These control samples 

are used to determine the correction factors which account for differences between the 

Pt 

data and simulation for high-p^^®® events. The factors are 1.03 ± 0.03 (1.38 ± 0.10) for the Th 
(light-lepton flavor) channels. This mismodeling arises from the mismeasurement of pj for jets 
and leptons. 

High-purity samples of Z —> tt —?■ events, from which the terms can be evaluated 

for the search channels with at least one Th, are obtained by removing the VBF selection and 
requiring mT(^, < ISGeV. The VBF efficiency for Z —> rr+jets processes is obtained 

from data using the Z —> £f'+jets control samples described above. 

The QCD multijet background is negligible for all final states except the ThThjj channel. To 
estimate the QCD multijet contribution to this channel, we select a QCD-dommated (>90% 
purity) CRl by requiring two Th candidates with relaxed isolation requirements. In addition, 
we require that the CRl events contain an SS ThTh pair. The SS signal region is thus included 
in CRl, but the potential impact of signal events is found to be negligible. The QCD multijet 
background in the OS ThTh channel is estimated by: 

'^QCD("^i7 ) = NgJ) (central) Rqs/ss ^vbf(^/7 ), (3) 

where (central) is the yield observed in the CRl sample with no VBF requirements, fol¬ 
lowing subtraction of the non-QCD component from CRl using simulation. The OS-to-SS ratio 
l^os/ss is obtained from a low-p™^®® (p^^®® < 30CeV) region after subtraction of the non-QCD 
contributions: we find Rqs/ss = 1-33 ± 0.03. Besides its use in the background determination 
procedure [Eq. Q], the measured result for Rqs/ss is used to provide a cross-check: we use it 
to extrapolate from an SS to an OS control region, both selected requiring p“'®® < 30 CeV and 
two non-isolated Th candidates. The obtained prediction for the rate of non-isolated OS Th lep¬ 
tons is in agreement with the observation. Finally, the efficiency is measured in exclusive 
sidebands fulfilling inverted Th isolation criteria. It is estimated as the rate of events with two 
non-isolated Th candidates plus two jets satisfying the VBF requirements divided by the rate 
of events with two non-isolated Th candidates without any additional jet requirements. The 
measured efficiency is evBF = 0.35% ± 0.08% (stat) ± 0.06% (syst). 

The QCD multijet background in the SS ThTh// channel is estimated using the following relation: 


. rpred _ yrfail-VBF 
^''QCD “ ^''QCD 


^non-isoTh 
-VBF 

-| non-isoTh' 
^ “ ^VBF 


(4) 


Here, is the observed yield in data, with non-QCD background from simulation sub¬ 

tracted, in an SS ThTh control sample requiring at least two jets not associated with one of the 
Th candidates to fail any of the |Ap|, r]ip 2 , or rrijj requirements. The VBF efficiency 
(cvBF for short) is measured in six exclusive Th isolation sidebands, without a p“'®® requirement 
and at least two jets. The validity of the method is demonstrated in data by the agreement that 
is observed, within statistical uncertainties, of these six independent measurements of Cvbf- 
The six corresponding control samples in simulation are used to test the stability of Cvbf as a 
function of pJp’®® and the Th isolation requirements. For this purpose, the probability for a sin¬ 
gle jet to be misidentified as a Th lepton is determined from simulation. The misidentification 
rates depend on the jet px and are used to determine an overall event weight by randomly 
selecting two jets in the event to represent the Th leptons. The VBF efficiencies in simulation 
are calculated from these weighted samples and demonstrate consistency with respect to the 
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7 Systematic uncertainties 


different and Th isolation requirements at the level of ~19%, which is assigned as a sys¬ 
tematic uncertainty in the background prediction. The VBF efficiency for nijj > 250 GeV is 
CVBF = 6.7% ± 0.5% (stat)lo; 5 o;;° (syst). 

7 Systematic uncertainties 

The main contributions to the total systematic uncertainty in the background predictions arise 
from the closure tests and from the statistical uncertainties associated with the data control re¬ 
gions used to determine the eyEF' (central), and Rqs/ss factors. The relative systematic 
uncertainties in SF^^^ and Rqs/ss related to the statistical precision in the CRs range between 
1 and 25%, depending on the background component and search channel. For nijj > 250 GeV, 
the statistical uncertainties in lie between 3 and 21%, while the systematic uncertainties 
evaluated from the closure tests and cross-checks with data range from 2 to 20%. For the back¬ 
ground eyfp, we assign no uncertainty due to the jet energy correction, as the mjj distributions 
are taken directly from the data control regions. 

Less significant contributions to the systematic uncertainties arise from contamination by non- 
targeted background sources to the GRs used to measure eyfp, and from the uncertainties in 
SFpQ^ (central) due to the lepton identification efficiency, lepton energy and momentum scales, 
pmi®® and trigger efficiency. 

The efficiencies for the electron and muon trigger, reconstruction, identification, and isolation 
requirements are measured with the "tag-and-probe" method 1^ ISTI with a resulting uncer¬ 
tainty of 2%. The Th trigger and identification-plus-isolation efficiencies are measured from a 
fit to the Z —> TT —> pTh visible mass distribution in a sample selected with a single-muon 
trigger, leading to a relative uncertainty of 4% and 6% per Th candidate, respectively Il46l . The 
p^^®® scale uncertainties contribute via the jet energy scale (2-10% depending on rj and px) 
and unclustered energy scale (10%) uncertainties, where "unclustered energy" refers to energy 
from a reconstructed object that is not assigned to a jet with px > 10 GeV or to a lepton with 
Px > 10 GeV. 

Since the estimate of the background is partly based on simulation, the signal and background 
rates are affected by similar sources of systematic uncertainty, such as the luminosity uncer¬ 
tainty of 2.6% HTI . The uncertainties in the lepton identification efficiency, lepton energy and 
momentum scale, pi^®® scale, and trigger efficiency also contribute to the systematic uncertainty 
in the signal. 

The signal event acceptance for the VBF selection depends on the reconstruction and identifi¬ 
cation efficiency and jet energy scale of forward jets. The jet reconstruction-plus-identification 
efficiency is >98% for the entire rj and px range, as is validated through the agreement ob¬ 
served between data and simulation in the t] distribution of jets, in particular at high rj, in 
control samples enriched with tf background events. The dominant uncertainty in the signal 
acceptance is due to the modelling of the kinematic properties of jets, and thus the VBF effi¬ 
ciency, for forward jets in the MadGraph simulation. This is investigated by comparing the 
predicted and measured nijj spectra in the Z+jets GRs. The level of agreement between the 
predicted and observed rrijj spectra is better than 23%, which is assigned as a systematic uncer¬ 
tainty in the VBF efficiency for signal samples. The uncertainty in the signal acceptance due to 
the PDF set included in the simulated samples is evaluated in accordance with the PDF4LHG 
recommendations HSj [13 by comparing the results obtained using the GTEQ6.6L, MSTW08, 
and NNPDFIO PDF sets 1141115011^ with those from the default PDF set (GTEQ6L1). The dom¬ 
inant uncertainties that contribute to the mjj shape variations include the pJp‘®®and jet energy 
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scale uncertainties. Correlations of the uncertainty sources are discussed in Section]^ 


8 Results and interpretation 

Figures|^and|^present the data in comparison to the predicted SM background. The combined 
results from all channels are shown in Fig.|^ Numerical results are given in Tables|^and|^ The 
observed numbers of events are seen to be consistent with the expected SM background in all 
search regions. Therefore the search does not reveal any evidence for new physics. 

Table 2: Number of observed events and corresponding background predictions for the OS 
channels. The uncertainties are statistical, including the statistical uncertainties from the con¬ 
trol regions and simulated event samples. 


Process 



a 

\ ^hii 

Z+jets 

4.3 ± 1.7 

3 7+^-^ 

a./_i 9 

19.9 ±2.9 

12.3 ± 4.4 

W+jets 

<0.1 

4 0+3.3 

17.3 ± 3.0 

2.0 ±1.7 

W 

2.8 ± 0.5 

3.1 ±0.7 

2.9 ±0.5 

0.5 ±0.2 

tt 

24.0 ± 1.7 

19.0^2.3 

11.7 ±2.8 

- 

QGD 

— 

— 

— 

6.3 ±1.8 

Higgs boson 

1.0 ±0.1 

1.1 ±0.5 

— 

1.1 ±0.1 

VBFZ 

— 

— 

— 

0.7 ±0.2 

Total 

32.2 ± 2.4 

o-i -I +4.6 

51.8 ±5.1 

22.9 ± 5.1 

Observed 

31 

22 

41 

31 


Table 3: Number of observed events and corresponding background predictions for the SS 
channels. The uncertainties are statistical, including the statistical uncertainties from the con¬ 
trol regions and simulated event samples. 


Process 


e±/i±;7 



Z+jets 

<0.1 

^-0 

0.5 ±0.2 

<0.1 

W+jets 

<0.1 

g +3.0 

^-0 

9.3 ±2.3 

0.5 ±0.1 

VV 

2.1 ±0.3 

1 q +0.4 

1.1 ±0.2 

0.1 ±6.5 X 10-2 

tf 

3.1 ±0.1 

3 5+‘^-^ 
a.a_o.9 

6.7 ±2.8 

0.1 ±1.2 X 10-2 

Single top 

— 


— 

<0.1 

QGD 

— 

— 

— 

7.6 ± 0.9 

Higgs boson 

— 

— 

— 

<0.1 

Total 

5.4 ± 0.3 

5 4±3-^ 

17.6 ±3.8 

8.4 ± 0.9 

Observed 

4 

5 

14 

9 


To quantify the sensitivity of this search, the results are interpreted in the context of the R- 
parity conserving minimal supersymmetric SM by considering production of charginos and 
neutralinos with two associated jets, as described in Section Models with a bino-like ^ 
and wino-like ^ and are considered. Since the ^ and Xi" belong to the same gauge group 
multiplet, we set ~ present results as a function of this common mass and the LSP 

mass . In the presence of a light slepton, I = e/ likely that the x^ will decay to £v^ 

and the ^ to The results are interpreted by considering £ = f and assuming branching 

fractions —)■ rf —> tt^) = 1 and B{x^ —)■ vt —)■ vrj^) = 1. To highlight how the VBF 

searches described in this paper complement other searches for the electroweak production of 
SUSY particles IITtIITsII . two scenarios are considered: (i) = OGeV (uncompressed-mass 

spectrum) and (ii) m~± — = SOGeV (compressed-mass spectrum). 
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Figure 3: Dijet invariant mass distributions in the (upper left) OS jiji, (upper right) SS 
(lower left) OS eji, and (lower right) SS eji signal regions. The signal scenario with = 
m^± = 200 GeV, rrif = 195 GeV, and = OGeV, as described in Section is shown. The 
signal events are scaled up by a factor of 5 for purposes of visibility. The shaded band in the 
ratio plot includes the systematic and statistical uncertainties in the background prediction. 

The cumulative signal event acceptance is shown in Table at three stages of the analysis: 
accounting for the branching fractions for the SUSY event to yield the indicated two-lepton 
channel (BF), the acceptance following application of the central selection (Gentral), and the 
acceptance following the VBF selection (VBF). The average px values of the e, p, and Th ob¬ 
jects in signal events are relatively soft, because of the energy and momentum carried by the 
associated neutrinos in the t decays. The OS and SS channels have similar signal acceptance 
because lepton pairs satisfying the event selection do not necessarily originate from the ^ or 
decays. The best signal sensitivity comes from the SS pp and ep channels due to a better 
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Figure 4: Dijet invariant mass distributions in the (upper left) OS (upper right) SS 
(lower left) OS and (lower right) SS signal regions. The signal scenario with = 

m^± = 200 GeV, nif = 195 GeV, and = OGeV, as described in Section is shown. The 
signal events are scaled up by a factor of 5 for purposes of visibility. The shaded band in the 
ratio plot includes the systematic and statistical uncertainties in the background prediction. 

background suppression with respect to a given signal acceptance. 

The expected signal yields from simulation with = OGeV and = 50GeV, are 

presented in Table The signal acceptance depends on the mass rrif of the intermediate r 
slepton. The results in Table are presented under two different assumptions for nif: (i) a 
fixed-mass difference assumption = 5GeV, and (ii) an average-mass assumption 

nif = 0.5m~± -|- 0.5?«^. In the compressed-mass-spectrum scenario, for which m~± — = 

50 GeV, the average-mass assumption yields significantly lower average lepton pj than the 






















































14 


8 Results and interpretation 


I 160; 

S 140- 

CM 


CMS 




500 


19.7 fb WS TeV) 


en, HT 


Others 

Signal (pp ^ xx\i) x 5 

(m_ =m-. = 200 GeV, m_ = 195 GeV, m_. = 0 GeV) 



1000 


1500 


2000 2500 

m(j,i) [GeV] 


Figure 5: Dijet invariant mass distribution for the combination of all search channels. The 
signal scenario with = 200 GeV, nif = 195 GeV, and = OGeV, as described in 

Section is shown. The signal events are scaled up by a factor of 5 for purposes of visibility. 
The shaded band in the ratio plot includes the systematic and statistical uncertainties in the 
background prediction. 

fixed-mass assumption, and the acceptance is lower by a factor of 2-3. In the uncompressed- 
mass-spectrum scenario, with ?«^=0GeV, the average-mass assumption produces larger av¬ 
erage lepton pj than the fixed-mass assumption, yielding an event acceptance that is 1.3-1.8 
times larger. 

Table 4: Gumulative signal event acceptance after application of the BF, central, and VBF re¬ 
quirements (see text). Note that the jet px threshold for the pp;; and ThTh;; final states is 30 GeV, 
while it is 50 GeV for the other final states. 


Gharmel 

BF(>l£i & >14) 

Gentral 

VBF 

p±T^^ (p±T±) 

0.399 

0.020 (0.020) 

0.007 (0.007) 

e^p^ (e^p^) 

0.152 

0.037 (0.037) 

0.014 (0.014) 


0.717 

0.010 (0.010) 

0.009 (0.009) 

(p^p^) 

0.081 

0.018 (0.018) 

0.007 (0.017) 


The calculation of the exclusion limit is obtained by using the mjj distribution in each channel 
to construct a combined likelihood in bins of nijj and computing a 95% confidence level (GL) 
upper limit on the signal cross section using the asymptotic GLg criterion Il52ti54l . Systematic 
uncertainties are taken into account as nuisance parameters, which are removed by marginal¬ 
ization, assuming a gamma or log-normal prior for normalization parameters, and Gaussian 
priors for mass spectrum shape uncertainties. The combination of the eight search channels 
requires simultaneous analysis of the data from the individual charmels, accounting for all sta¬ 
tistical and systematic uncertainties and their correlations. Gorrelations among backgrounds, 
both within a charmel and across channels, are taken into consideration in the limit calculation. 
For example, the uncertainties in physics object identification and reconstruction are treated as 
correlated for charmels with a common particle in their final states, while the uncertainty in 
the integrated luminosity is treated as correlated across channels. The uncertainties resulting 
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Table 5: Signal event yields from simulation. The first terms {m^±,mf} correspond to the fixed- 
mass difference assumption = 5GeV, while the terms in parentheses 

correspond to the average-mass assumption nif = 0.5m^ + 0.5m^±. 


{m{xf),m{f)} [GeV] 

(loose) 

(tight) 

mi 

Fmi 

'GTh;; 


m 

(^) = 0GeV 




{100, 95} ({100,50}) 

16(29) 

6.6(12) 

13(24) 

7.1 (9.4) 

8.7(10.7) 

{200,195} ({200,100}) 

5.4(9.7) 

1.8(3.1) 

3.5(6.3) 

4.5(6.0) 

3.8(4.7) 

{300, 295} ({300,150}) 

2.3(4.1) 

0.68(1.2) 

1.4(2.4) 

1.9(2.5) 

1.5(2.0) 

{400, 395} ({400, 200}) 

0.57(1.0) 

0.17(0.30) 

0.35(0.62) 

0.46(0.63) 

0.38(0.51) 

Am(^-^) =50GeV 

{200,195} ({200,175}) 

1.4(0.5) 

0.85(0.33) 

1.7(0.65) 

0.99(0.35) 

0.46(0.09) 

{300, 295} ({300, 275}) 

0.47(0.18) 

0.28(0.11) 

0.58(0.23) 

0.40(0.14) 

0.20(0.04) 

{400, 395} ({400, 375}) 

0.12(0.05) 

0.08(0.03) 

0.15(0.06) 

0.10(0.03) 

0.05(0.01) 


from the number of simulated events, and from the event acceptance variation with different 
sets of PDFs in a given mjj bin, are treated as uncorrelated within a channel and correlated 
across channels. The uncertainties due to the closure tests are treated as uncorrelated within 
and across the different final states. 

Figures (left) and (right) show the expected and observed limits as well as the theoretical 
cross section as functions of m~± for, respectively, the fixed- and average-mass mt assumptions. 
For the fixed-mass assumption with a compressed-mass spectrum {m^± — = SOGeV),^ 

and gauginos with masses below 170 GeV are excluded, where the previous ATLAS and 
GMS SUSY searches do not probe. For the average-mass assumption with an uncompressed- 
mass spectrum (iw^ = 0), the corresponding limit is 300 GeV. These mass limits are conserva¬ 
tively determined using the theoretical cross section minus its one standard deviation uncer¬ 
tainty. The limits of 320 and 380 GeV for = OGeV in Refs. llZlIIlI can be compared 
to the corresponding result of 300 GeV in the present analysis [see the yellow band in Fig. 
(right)]. 


9 Summary 

A search is presented for non-coloured supersymmetric particles in the vector-boson fusion 
(VBF) topology using data corresponding to an integrated luminosity of 19.7fb^^ collected 
with the GMS detector in proton-proton collisions at v^s = 8 TeV. This is the first search for 
SUSY in the VBF topology. The search utilizes events in eight different final states covering 
both same- and opposite-sign dilepton pairs. The leptons considered are electrons, muons, 
and hadronically decaying t leptons. The VBF topology requires two well-separated jets that 
appear in opposite hemispheres, with large invariant mass mjj. The observed mjj distributions 
do not reveal any evidence for new physics. The results are used to exclude a range of xf and 
^ gaugino masses. For models in which the ^ lightest-supersymmetric-particle mass is zero, 
and in which the xf and ^ branching fractions to T leptons are large, xf and ^ masses up to 
300 GeV are excluded at 95% GL. For a compressed-mass-spectrum scenario, in which m~± — 
m^^ = 50 GeV, the corresponding limit is 170 GeV. While many previous studies at the LHG 
have focused on strongly coupled supersymmetric particles, including searches for charginos 
and neutralinos produced in gluino or squark decay chains, and a number of studies have 
presented limits on the Drell-Yan production of charginos and neutralinos, this analysis obtains 
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9 Summary 
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Figure 6: Combined 95% CL upper limits on the cross section as a function of = m~±. The 
signal cross section is calculated with the VBF jet selection: jet px > 30 GeV, |Ap(jets)| > 4.2, 
and P 1 P 2 < 0. (left) The results for the fixed-mass difference assumption, in which m~± — nif = 
5 GeV, for m~± —m^ =50 GeV (compressed-mass spectrum) and m^= 0 GeV (uncompressed- 
mass spectrum), (right) The corresponding results for the average-mass assumption, in which 


nif = 0.5m 


xt 


■ 0.5m^. 

Xi 


the most stringent limits to date on the production of charginos and neutralinos decaying to r 
leptons in compressed-mass-spectrum scenarios defined by the mass separation Am = m^± — 


m^ < 50 GeV. 
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